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SUMMARY 

In  visual  psychophysics  it  is  sometimes  necessary  to  be  able  to  display  lumi¬ 
nances  with  12  bit  precision.  This  report  describes  a  method  where  two  8  bit 
color  outputs  of  a  DeAnza  IP8400  image  processor  are  combined  into  a  mono¬ 
chrome  signal  with  an  effective  resolution  of  12  bits.  The  implementation  in 
hardware  and  software,  and  the  calibration  of  the  system  are  described.  Measur¬ 
ements  show  that  the  system  meets  the  desired  specifications. 
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Rap. nr.  IZF  1992  B-7 


Instituut  voor  Zintuigfysiologie  TNO 

Soesterberg 


Hoge-resolutie  contrastweergave  op  een  video  monitor:  methode  en  ijking 
J.M.  Valeton  en  A.J.C.  de  Reus 

SAMENVATTING 

Bij  visuele  psychofysica  is  soms  een  luminantie-resolutie  van  12  bit  noodzakelijk. 
In  dit  rapport  wordt  een  methode  uitgewerk’.,  waarbij  twee  8  bit  kleurenuit- 
gangen  van  een  DeAnza  IP8400  image  processor  gecombineerd  worden  tot  een 
monochroom  signaal  met  een  effectieve  resolutie  van  12  bit.  De  implementatie 
in  hardware  en  software  wordt  behandeld,  alsmede  de  ijking  van  het  systeem. 
Uit  metingen  blijkt  dat  het  systeem  aan  de  verwachte  prestaties  voldoet. 
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1  INTRODUCTION 

In  experiments  on  contrast  detection  it  is  often  necessary  to  present  spatially 
complex  spatial  stimuli  at  very  low  contrast  (as  low  as  0.5%),  with  a  specific 
temporal  profile.  Displaying  such  stimuli  accurately  on  a  display  system  that  has 
basically  an  8  bit  intensity  resolution  is  not  possible.  In  this  report  we  will  first 
show  that  an  effective  resolution  of  at  least  12  bit  is  required.  In  the  rest  of 
Chapter  2  we  describe  how  a  method,  proposed  by  Watson  et  al.  (1986)  to  solve 
this  problem,  was  implemented  on  a  DeAnza  IP8400  image  processor  with  a 
Gould  Concept  32/67  host  computer.  Two  (8  bit)  channels  of  the  RGB  color 
output  of  the  image  processor  are  combined  to  drive  a  monochrome  display,  and 
hence  the  method  is  called  two  channel  contrast  control.  In  Chapter  3  we  will 
describe  how  such  a  system  is  calibrated.  In  Chapter  4  we  show  the  results  of 
test  measurements  to  check  the  performance  of  the  two  channel  method.  In  the 
appendices  the  source  code  of  our  12  bit  display  routine  and  the  technical  data 
of  the  video  monitor  used  in  the  test  measurements  are  presented. 

In  some  parts  the  text  of  this  report  relies  heavily  on  the  paper  by  Watson  et  al. 
(1986).  The  report  forms  the  theoretical  background  of  12  bit  stimulus  display  in 
our  laboratory.  Our  implementation  of  Watson’s  method  and  the  calibration 
procedures  are  documented  in  a  separate  report  as  part  of  the  laboratory  quality 
control  procedures  (De  Reus  &  Valeton,  1992). 

In  the  report,  parameters  that  appear  in  the  computer  programs  are  printed 
within  square  brackets,  e.g.  [red  offset].  Software  algorithms  are  presented  in 
pseudo-code. 


2  TWO  CHANNEL  CONTRAST  CONTROL 
2.1  The  12-bit  resolution  requirement 

A  digital  display  system  renders  images  in  discrete  intensity  steps  and  hence  is 
prone  to  quantization  errors.  For  accurate  rendering  of  a  sinewave  a  minimum 
of  8  to  10  intensity  levels  is  required.  A  display  system  with  a  b  bit  digital-to- 
analog  converter  (DAC)  has  a  smallest  voltage  step  size1  of  l/(2b  -  1).  The 
corresponding  lowest  possible  (voltage)  contrast  of  a  sinewave  on  a  constant 
background  is  2/(2b  -  /),  which  for  an  8  bit  system  amounts  to  0.78%.  At  this 
contrast,  however,  there  are  only  2  intensity  levels  available,  so  the  signal  is  a 
square  wave  instead  of  a  sinewave.  The  contrast  threshold  for  spatial  sinewave 
patterns  is,  under  optimal  conditions,  on  the  order  of  0.5%.  In  threshold  experi¬ 
ments.  therefore,  sine  waves  must  be  reproduced  accurately  at  contrasts  well 


A  system  with  a  b  bit  output-DAC  has  2h  signal  levels.  Thus  an  8  bit  system  has  256  different 
levels. 
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below  0.5%.  Since  at  0.78%  contrast  an  8  bit  system  is  already  limited  to  only  2 
levels,  it  is  clear  that  an  8  bit  system  lacks  sufficient  intensity  resolution.  The 
problem  is  illustrated  in  Fig.  1. 


Fig.  1  Sinewave  with  0.3%  contrast  rendered  at  different  uniform 
quantization  levels.  The  solid  line  is  the  original  signal,  the  dashed 
line  is  the  sinewave  quantized  to  256  levels  (8  bits),  the  broken  line  to 
776  levels  (9.6  bits)  and  the  dotted  line  to  4096  levels  (12  bits).  Only 
at  full  contrast  (100%)  all  levels  are  used. 

In  this  figure,  taken  from  Pelli  and  Zhang  (1991),  we  see  what  happens  when  a 
0.3%  contrast  sinewave  (solid  line)  is  displayed  with  an  8  bit,  a  9.6  bit  and  a  12 
bit  resolution  system.  An  8  bit  system  results  in  zero  output,  an  effective  9.6  bit 
output  generates  3  intensity  levels  and  a  12  bit  DAC  gives  a  reasonable  electron¬ 
ic  rendering  of  the  sinewave.  This  shows  that  for  sinewave  like  signals  at 
threshold  contrast  levels,  a  12  bit  resolution,  or,  a  luminance  step  size  of  about 
1/4000  (of  full  scale)  is  required. 

Such  a  high  intensity  resolution  output  for  a  display  can  theoretically  be  obtained 
by  using  12  bit  DACs,  but  such  a  system  is  not  on  the  market  and  building  one 
would  be  very  costly.  An  elegant  and  simple  alternative,  called  "nvo  channel 
contrast  control",  was  proposed  by  Watson  et  al.  (1986).  Pelli  and  Zhang  (1991) 
presented  an  extension  of  this  approach  to  three  channels. 

In  this  method,  two  channels  of  the  8  bit  RGB  color  output  of  a  display  system 
are  combined  to  drive  a  monochrome  monitor.  One  channel  is  attenuated  by  a 
certain  factor  and  added  to  the  other  channel.  The  attenuated  channel  can 
display  only  part  of  the  whole  intensity  range  but  it  has  all  8  bits  available.  This 
results  in  a  smaller  step  size  and  hence  a  higher  intensity  resolution.  The  other 
channel  is  added  to  put  the  high  resolution  section  on  the  required  average 
intensity  level.  For  large  contrasts  only  the  unattenuated  channel  is  used  and 
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therefore  in  this  case  the  step  size  is  not  smaller  than  with  a  normal  8  bit  display 
system. 


2.2  Hardware  implementation 

The  output  voltage  of  the  red  and  green  channels,  VR  and  VG,  of  the  DeAnza 
image  processor  are  added  electronically  to  yield  the  following  output  voltage  V: 

V  -  (1  -f)  VR  +  f  vG 

where  /  is  an  attenuation  factor  with  0  <  f  <  1.  We  chose  /  =  1/16.  The  setup  is 
illustrated  in  Fig.  2. 


voc 


I  l 

R-lut 


Fig.  2  Hardware  diagram  of  the  video  output  controller  (VOC)  of  the 
DeAnza  image  processor  and  the  two  channel  addition. 


In  the  green  channel  the  step  size  is  decreased  by  a  factor  of  ///,  which  means 
that  the  intensity  resolution  is  increased  bv  the  same  factor  (16.  or:  4  bits).  The 
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total  resolution  of  the  combined  system,  for  the  limited  intensity  range,  now 
effectively  becomes  12  bits. 

The  Video  Output  Controller  (VOC)  is  the  hardware  through  which  images  in 
memory  are  sent  to  the  CRT.  There  are  DACs  for  each  color  channel  R,  G,  and 
B  and  there  is  also  a  channel  for  Graphic  Overlay.  A  Look-Up  Table  (LUT) 
precedes  each  DAC,  and  each  LUT  has  four  sections,  numbered  0  -  3.  These 
sections  each  have  256  entries  and  are  fully  equivalent.  A  LUT  section  can  be 
assigned  by  software.  Different  signal  paths  in  the  VOC  can  be  set  up  under 
software  control.  In  the  present  situation  the  memory  plane  containing  the 
stimulus  image  is  connected  to  the  R,  G,  and  B  DACs,  through  LUT  sections  0; 
the  Graphic  DAC  is  not  used.  The  B-LUT  (blue  channel)  is  set  to  zt.o,  so  this 
channel  is  effectively  not  connected.  The  memory  planes  in  the  DeAnza  each 
have  an  8  bit  LUT  of  their  own  (called  Intensity  Transformation  Table,  or  ITT) 
but  these  are  bypassed  in  the  present  setup. 

The  DeAnza  IP8400  display  system  has  8  bit  resolution  output.  As  Pelli  and 
Zhang  (1991)  point  out,  8  bit  DACs  have  an  effective  resolution  of  only  7  bits 
due  to  noise  in  the  least  significant  bit.  In  the  DeAnza  this  problem  has  been 
elegantly  solved  by  building  10  bit  DACs  into  the  VOC  and  then  using  the  8 
most  significant  bits  only. 


2.3  Software  implementation 

The  maximum  contrast  that  can  be  displayed  through  the  attenuated  channel  is 
called  the  "critical  contrast"  Ccnr  For  a  linear  monitor  the  critical  contrast  would 
be  equal  to  /  (1/16  in  our  implementation),  but  for  practical  CRT  displays  the 
critical  contrast  is  slightly  larger  than  /  (up  to  a  factor  of  2,  see  Chapter  3).  The 
system  operates  in  two  modes,  for  contrasts  below  and  above  Ccri(  respectively. 

Mode  1  Contrasts  higher  than  Ccri[.  Only  the  R  channel  is  used  and  the 

G-LUT  is  filled  with  zeros.  This  mode  is  in  fact  a  traditional 
single  channel  system.  Contrast  is  manipulated  by  changing  the 
R-LUT  (red  channel). 

Mode  2  Contrast  lower  than  Ccnt.  Two  channels  are  used.  Signal  contrast 

is  manipulated  by  changing  the  G-LUT  (green  channel)  while 
the  R-LUT  is  filled  with  the  constant  [red_offset].  This  mode 
generates  a  limited  range  of  contrasts,  but  has  an  effective 
intensity  resolution  of  12  bit. 

The  value  [red_offset]  is  chosen  such  that  when  zero  contrast  must  be  displayed, 
and  the  G-LUT  is  filled  with  [bkg_value]  (=  e.g.  128),  the  monitor  luminance  is 
equal  to  the  lumincnce  in  mode  1  when  [bkg  value]  is  loaded  into  the  R-LUT 
(to  produce  the  stimulus  background).  The  following  equation  roughly  holds: 
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L2  {  R-LUT[red_offset]  +  G-LUT[bkg_value]  )  =  L,  {  R-LUT[bkg_value]  } 

where:  L[  and  U,  are  the  screen  luminances  in  mode  1  and  2.  The  procedure  for 
determining  [red_offset]  is  part  of  the  system  calibration  and  is  given  in  detail  in 
Chapter  3.  When  [red_offset]  is  chosen  carefully,  the  mean  display  and  back¬ 
ground  luminances  do  not  change  when  the  system  switches  between  mode  1  and 
mode  2. 

The  source  code  for  the  routine  which  implements  the  12  bit  display  using  two 
modes  is  given  in  Appendix  1.  Within  the  subroutine  CMVARF2()  there  are  two 
main  streams.  When  the  requested  contrast  is  greater  than  Ccrit,  a  section  of 
code  according  to  mode  1  is  executed,  and  when  the  contrast  is  smaller  than 
Ccrit,  the  two  channel  system  of  mode  2  is  used.  Note  that  in  the  latter  case  the 
contrast  is  divided  by  Ccrit.  This  is  where  the  increase  in  contrast  resolution 
actually  happens:  by  multiplying  a  small  contrast  with  1/Ccrit,  there  are  about  1/f 
times  more  levels  available  for  rendering  the  signal.  At  the  output  the  signal  is 
multiplied  by  f  in  analogue  hardware.  The  result  thus  is  more  available  lumi¬ 
nance  levels,  but  only  in  a  small  luminance  range. 


2.4  Dynamic  contrast  control 

In  this  section  we  explain  how  the  contrast  of  the  stimuli  on  the  display  is 
controlled.  Correction  for  the  non-linearity  of  the  CRT  display  is  treated  in  the 
next  section  and  in  Chapter  3. 

The  image  that  must  be  displayed  is  stored  in  an  8  bit  512x512  pixel  memory 
channel  at  maximum  contrast,  i.e.  all  256  available  levels  are  used.  This  is  done 
to  have  the  maximum  number  of  levels  available  in  the  output  signal  in  order  to 
minimize  the  quantization  error  in  the  image  itself.  In  many  cases  the  signal  to 
be  displayed  is  smaller  than  512x512  pixels.  In  those  cases  the  signal,  e.g.  a  patch 
of  a  sinewave  grating,  is  usually  placed  in  the  center  of  the  image  and  the 
surrounding  image  area  is  set  to  0.  We  reserve  this  pixelvalue  to  display  the 
background  intensity  using  the  LUT.  The  8  bit  LUTs  each  have  256  entries.  The 
first  position  of  the  look-up  table  is  used  to  set  the  background  that  surrounds 
the  signal,  or:  LUT[0]  =  [bkg_index],  where  [bkg_index]  is  usually  set  to  128.  An 
uneven  number  of  255  LUT  positions  remain  for  use  in  setting  the  signal 
contrast.  The  middle  position  is  128  and  there  are  127  positions  above  and  below 
this  level.  This  enables  us  to  generate  both  normal  and  reversed  contrasts 
symmetrically  around  zero. 

The  contrast  of  the  signal  on  the  display  is  manipulated  by  "rotating"  the  LUT 
contents  around  the  center  position  S  ( =  128),  see  Fig.  3. 
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Fig.  3  Illustration  of  the  rotating  LUT.  Line  1  (dotted)  represents  the 
LUT  transformation  for  a  full-contrast  (100%)  stimulus  display,  line  2 
(dashed)  the  LUT  transformation  for  a  negative  100%  contrast.  The 
tangent  of  the  angles  gives  contrasts  between  0  and  100%,  signed 
positive  for  a,  and  negative  for  a2- 
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A  LUT  loaded  with  a  linear  ramp  (solid  line)  produces  an  image  of  maximum 
contrast,  a  LUT  with  the  value  128  at  each  location  yields  a  zero  contrast  "grey" 
image  and  a  LUT  with  a  slope  of  -1  (dashed  line)  shows  the  image  with  reversed 
contrast.  Rotation  of  the  LUT  around  point  S  effectively  changes  the  contrast  to 
any  desired  level:  the  contrast  is  equal  to  the  slope  (or  tan(a),  where  -t/4  <  a 
<  t/4  when  S  =  128)  of  the  line  in  Fig.  3. 

A  stimulus  image  can  be  displayed  dynamically,  with  a  specified  time  course,  by 
loading  a  new  LUT  for  every  video  frame.  In  practice  this  is  implemented  by 
using  a  list  of  numbers  [tjist],  e.g.  a  Gaussian  profile  as  in  Fig.  4.  representing 
the  contrast  for  each  frame. 
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Fig.  4  Gaussian  time  course  for  dynamic  contrast  display.  This 
example  extends  1  second  at  66  fields/second. 


t  jist  (gaussian) 


T  .  ■  .T-r,  .  .  ■  ■  n%ij 


To  create  the  contrasts  for  a  Gaussian  profile,  the  LUT  is  rotated  around  point 
S  in  Fig.  3  from  a  horizontal  line  (slope  0)  up  to  slope  1  and  back  down  again. 
For  maximum  resolution  the  time  profile  is  stored  at  full  contrast,  so  in  order  to 
display  the  stimulus  at  a  contrast  smaller  than  1.0  the  slope  of  the  LUT  for  each 
frame  must  be  multiplied  by  the  desired  contrast  (see  the  source  code  of  the 
display  routine  in  Appendix  1).  For  this  type  of  dynamic  stimulus  display  the 
system  must  be  fast  enough  to  recalculate  and  load  several  LUTs  between  two 
video  frames.  The  Gould/DeAnza  is  able  to  calculate  and  fill  two  LUTs  between 
two  fields  at  a  display  rate  of  66  fields/second.  (In  the  mode  we  use.  the 
DeAnzas  line-scanning  system  is  interlaced:  each  frame  is  composed  of  two 
fields.) 


2.5  Display  linearization 

As  will  be  explained  fully  in  Chapter  3,  a  provision  must  be  made  to  correct  for 
the  nonlinearity  of  the  CRT  display.  This  can  be  done  bv  using  a  correction  table 
as  look-up  table.  However,  the  look-up  table  is  also  used  for  controlling  the 
stimulus  contrast.  In  joftware  both  these  function  are  carried  out  at  the  same 
time,  using  a  single  look-up  table,  as  follows. 

For  each  image  pixelvaiue  i,  given  the  time  course  and  desired  contrast,  the 
value  is  calculated  as  if  the  monitor  were  linear,  using  the  mechanism  of  rotating 
a  LUT.  In  order  to  correct  for  the  display  non-linearity  these  new  values  are 
used  as  indices  for  the  array  correction  table  CORLUT  for  filling  the  actual 
hardware  display  LUT: 

VOC_LUT[i]  =  CORLUT[  ROTATE[i]  ] 
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2.6  Outline  of  the  subroutine  CMVARF2 

The  subroutine  CMVAR2F  takes  care  of  displaying  stimuli  of  a  specified  time 
course  using  the  two  channel  method.  It  is  written  in  Fortran  because  on  the 
Gould  computer  we  use  Fortran  is  proven  to  be  faster  than  C.  Additional 
advantages  are  that  the  low-level  DeAnza  interface  is  written  in  Fortran  and  that 
most  of  the  software  used  on  the  Gould/De Anza  system  is  in  Fortran.  Fortran 
subroutines  and  functions  can  be  called  from  C  programs  but  not  vice  versa. 

Subroutine  header: 

SUBROUTINE  CMVARF2(USER,  TLIST,  ATTEN,  B  KG  INDEX,  CHAN,  CO, 
NFRAMES,  MODE,  CRIT  CONT,  CCM1,  CCM2) 


where: 


Parameter 

Fortran  tvpe 

C  type 

Description 

USER 

INTEGERM 

int 

DeAnza  user  number  for  multi-user  environ¬ 
ment 

TLIST 

REALM  ARRAY 

float  * 

time  course,  first  element  is  length  of  this 
table 

ATTEN 

REALM 

float 

stimulus  attenuation 

BKGINDEX 

INTEGERM 

int 

background  index  in  CCM1,  used  for  all  pixels 
with  value  0 

CHAN 

INTEGERM 

int 

memory  plane  number  of  stimulus 

CO 

INTEGERM 

int 

- 

NFRAMES 

INTEGERM 

int 

number  of  frames  in  the  stimulus 

MODE 

INTEGERM 

int 

display  mode:  1  =  8  bit,  2  =  12  bit 

CRIT  CONT 

REALM 

float 

critical  contrast  for  12  bit  mode 

CCM1 

INTEGERM  ARRAY 

int  * 

gamma  correction  table  for  red  channel 

CCM2 

INTEGERM  ARRAY 

int  * 

gamma  correction  table  for  green  channel 

The  Fortran 

and  C  data  types  are 

for  the 

compilers  installed  on  our  Gould  CONCEPT/32 

machine.  The  parameter  CO  is  not  for  user-purposes. 


The  subroutine  is  split  in  two  major  parts,  one  for  the  8  bit  mode  and  one  for 
the  two  channel  12  bit  mode.  The  8  bit  part  has  no  specifically  interesting 
fea  ares.  In  the  12  bit  mode,  for  each  frame  the  contrast  is  determined  (from  the 
desired  maximum  contrast  and  the  current  value  in  the  time  course).  When  this 
contrast  is  below  the  critical  contrast  the  green  LUTs  are  rotated  and  the  red 
LUT  is  set  at  [red  offset],  in  the  other  case  the  green  channel  is  set  to  zero  and 
the  red  LUT  is  rotated.  The  subroutine  is  optimized  for  speed  and  two  LUTs 
can  jc  calculated  and  filled  within  one  field  time  (1/66  second).  The  complete 
source  code  is  presented  in  Appendix  1. 
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3  SYSTEM  CALIBRATION 
3.1  CRT  non-linearity 

A  CRT  monitor  has  a  basically  non-linear  relation  between  driving  signal  v  and 
luminance  output  L  that  can  roughly  be  described  by  a  power  function: 

L  =  L  +  a  vY  (1) 

mui 

where  Lrajn  is  the  luminance  at  zero  input  and  gamma  usually  between  2  and  3. 
This  function  is  often  called  the  "gamma  function",  after  the  exponent  y.  The 
effect  of  this  nonlinearity  is  that  an  image  stored  in  the  display  memory  is  not 
rendered  truthfully  on  the  display.  For  most  psychophysical  stimuli  a  linear 
relation  between  pixel  value  and  display  luminance  is  required  and  this  can  be 
realized  by  displaying  the  stimuli  through  a  lookup  table  (LUT)  that  contains  the 
inverse  of  L(v).  When  "two  channel  contrast  control"  is  used,  the  situation  is 
slightly  more  complex  because  two  LUTs  have  to  be  calculated  and  loaded.  The 
calibration  program  Calibra  collects  the  L(v)  curve  for  two  channels,  the 
calculation  of  the  correction  LUTs  and  the  determination  of  [red  offset]  and 
Cent'  tw0  parameters  that  are  used  in  the  two  channel  method.  In  the  following 
we  will  present  step-by-step  procedures  to  determine  the  various  quantities. 


3.2  Measuring  the  gamma  function 

The  gamma  function  is  measured  in  an  automated  fashion  by  measuring  the 
display  luminance  with  a  Pritchard  Photometer  for  all  256  pixel  values.  The 
measurement  procedure  is  computer  controlled. 

3.2.1  Choice  of  test  pattern 

In  most  monitors  the  high  tension  anode  voltage  (e.g.  ca.  25  kV)  is  not  stabilized 
and  as  a  consequence  this  voltage  drops  with  increasing  beam  current.  This 
means  that  the  luminance  generated  by  a  patch  of  certain  (large)  pixel  value,  e.g. 
200,  is  smaller  for  a  large  patch  than  for  a  small  patch.  This  would  make  a 
measurement  of  the  gamma  function  dependent  on  stimulus  area  and  this  is 
obviously  undesirable.  To  more  or  less  circumvent  this  problem  we  use  test 
images  for  the  calibration  that  consist  of  2  areas,  see  Fig.  5. 
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Fig.  5  Three  calibration  patterns.  Each  pattern  contains  5122  pixels. 
The  white  areas  in  the  patterns  contain  pixels  with  value  0,  the  black 
areas  pixels  with  value  255.  The  gray  part  in  the  rightmost  pattern 
contains  pixels  with  value  128. 


One  area  of  the  first  two  test  patterns  is  light  (pixel  value  255)  and  the  other 
area  is  dark  (pixel  value  0).  Both  these  areas  are  varied  during  the  measure¬ 
ments  with  Calibra.  The  areas  (number  of  pixels)  of  the  light  and  dark  regions 
are  equal,  both  for  the  bar  pattern  and  for  the  square  pattern.  If,  when  the  pixel 
value  of  the  dark  area  is  increased,  the  pixel  value  for  the  light  area  is  decreased 
at  the  same  rate,  the  load  on  the  CRT  is  constant,  and  the  dependence  of  lumi¬ 
nance  on  pattern  size  is  removed.  In  a  third  test  pattern  we  replaced  the  light 
area  with  a  border  of  a  constant  grey  (pixel  value  128).  This  does  not  satisfy  the 
demand  of  a  constant  load  and  can  therefore  only  be  used  for  CRTs  with  a 
sufficiently  stabilized  anode  voltage.  The  bar  pattern  has  the  additional  advan¬ 
tage  over  the  square  patterns  that  the  average  load  for  each  scanline  is  constant. 

A  second  factor  in  the  choice  of  test  pattern  is  stray  light.  Light  is  scattered 
within  the  phosphor  and  in  the  glass  monitor  face.  When  a  dark  bar  is  flanked 
on  both  sides  by  a  bright  bar,  stray  light  from  the  bright  bars  precludes  a  correct 
luminance  measurement  for  the  dark  bar.  For  this  reason  the  dark  area  should 
be  as  large  as  possible.  The  square  test  pattern  has  the  advantage  that  the  bright 
surround  is  removed  as  far  as  possible  from  the  center  of  the  dark  area.  Other 
problems  with  stray  light  can  arise  when  light  not  coming  from  the  measurement 
part  of  the  test  pattern  falls  on  the  measurement  spot  of  the  photometer.  Care 
should  be  taken  to  avoid  problems  of  this  kind. 

3.2.2  Procedure  to  record  the  gamma  function 

The  procedure  to  record  the  gamma  function  is  as  follows: 

1  Load  a  light-dark  test  pattern  into  display  memory  (dark  area  has  pixel  value 
0;  light  area  has  pixel  value  255). 

2  Aim  the  photometer  at  the  center  of  the  dark  area. 

3  Set  the  lookup  table  LUT[0]  =  0  and  LUT[255]  =  255. 

4  Repeat 

5  Record  the  luminance  value. 
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6  Change  the  lookup  table  entries  as  follows: 

-  LUT[0]  =  LUT[0]  +  1, 

-  LUT[255]  =  LUT[255]  -  1. 
until  LUT[0]  equals  255. 

Note  that  the  choice  of  the  LUT  indices  0  and  255  is  arbitrary. 

3.2.3  Signal  averaging 

Luminance  is  recorded  with  a  Spectra  Physics  Pritchard  photometer.  This 
instrument  is  equipped  with  an  TNO-IZF  custom  made  24  bit  digital  interface  (a 
so-called  Normalized  Output  Connector,  or  NOC).  This  interface  is  connected  to 
NIC  (a  Normalized  Input  Connector)  on  the  Gould  host  computer.  The 
Pritchard’s  NOC  sends  out  3  luminance  readings  per  second  to  the  host 
computer.  In  order  to  get  an  accurate  estimate  of  the  display  luminance  a 
number  of  readings  are  averaged  until  an  error  limit  or  a  maximum  number  of 
samples  is  reached.  The  following  procedure  is  used: 

1  Set  pixelvalue  to  desired  value  (by  setting  the  LUT). 

2  Start  reading  luminance  values  and  discard  first  few  values  (get  rid  of  errors 
due  to  range  switching  and  let  the  photometer  adapt  to  the  new  luminance). 

3  Initialize  running  average  (set  to  zero)  and  set  COUNT  to  0 

4  Repeat 

5  Read  luminance  value  from  photometer. 

6  Update  running  average,  measurement  variance  and  COUNT. 

until  (COUNT  >  MIN  NUM  and  variance  <  ERROR  LIMIT)  or  (COUNT 
>  MAX  NUM). 

7  If  (variance  >  ERROR  LIMIT) 

8  Give  a  prompt  to  re-adjust  the  setup. 

9  Restart  this  routine  again. 

10  Return  the  running  average  as  estimate  for  the  screen  luminance. 

Usually  10  (=  MIN_NUM)  luminance  values  are  recorded  for  one  estimate  and 
the  error  limit  is  set  to  1%  of  the  running  average.  A  complete  calibration  run, 
for  a  single  channel,  takes  about  10  minutes.  An  example  of  a  gamma  curve  is 
presented  in  Fig.  6  as  G(d). 
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gamma  function 


pixelvalue 


Fig.  6  Example  of  a  measured  gamma  function. 


3.3  Linearization  of  a  single  channel  display 
3.3.1  Notation 

In  the  subsequent  sections  we  use  the  following  notation: 
p  digital  image  pixel  value  (0  =  <  p  =  <  255);  is  also  a  LUT  address', 
pO  central  pixel  value  (128) 

d  digital  value  of  a  LUT  entry  (0  =  <  d  =  <  255);  sent  to  the  CRT 
dO  central  LUT  entry  (128) 

Lmax  measured  luminance  for  pixel  value  255 

Lmin  measured  luminance  for  pixel  value  0  (zero  input) 

L0  background  display  luminance  (stimuli  are  superimposed  on  this  value) 
G(d)  gamma  function  (as  a  function  of  LUT  value  d) 

L(p)  function  representing  a  linear  relation  between  image  pixel  value  (p)  and 
luminance  (1  <,  p  <  255) 

C(p)  correction  table,  inverted  gamma  function. 

In  Fig.  7  most  of  these  quantities  are  shown. 


Fig.  7  Graphical  demonstration  of  linearizing  a  CRT.  See  the  text  in 
sections  3.3.1  and  3.3.4  for  details. 


Fig.  7  illustrates  how  the  correction  table  C(p)  is  constructed.  The  upper  left 
quadrant  shows  the  gamma  function  G(d)  as  it  is  actually  measured:  recorded 
luminance  as  a  function  of  LUT  entry.  The  minimum  luminance  Lmin  (>  0 
cd/m2)  is  obtained  for  LUT  entry  0,  Lmax  for  LUT  entry  255.  The  desired 
relation  L(p)  is  shown  in  the  upper  right  quadrant:  a  linear  relation  between 
image  pixel  values  and  screen  luminance,  (note  that  a  different  desired  relation 
can  also  be  used.)  Since  p  =  0  is  used  to  produce  the  stimulus  background 
luminance  L0,  L(l)  (and  not  L( 0 ) ! )  equals  Lmin;  L(255)  equals  Lmax.  The  image 
pixel  values  are  identical  to  LUT  addresses,  which  is  shown  by  the  mirror 
operation  in  the  lower  right  quadrant.  Finally,  the  lower  left  quadrant  shows  how 
the  output  LUT  behaves:  it  uses  image  pixel  values  p  as  addresses  and  translates 
these  p’s  into  d’s,  the  LUT  entries.  These  d’s  are  converted  to  analog  values 
(voltages)  by  the  output  DAC  and  are  displayed  on  the  screen. 

3.3.2  Procedure  to  determine  the  background  luminance  L0 

In  pseudo-code: 

1  Set  green  LUT  to  0. 

2  Load  the  value  of  255  into  the  red  LUT. 

3  Record  the  screen  luminance  as  Lmax. 

4  Load  the  value  of  0  into  the  red  LUT. 

5  Record  the  screen  luminance  as  Lmin. 

6  Compute  the  mean  luminance  LO  =  (Lmax  +  Lmin)  /  2. 
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This  procedure,  however,  is  not  used  in  the  present  setup.  We  calculate  L0  from 
the  gamma  function  of  the  red  channel,  according  to  step  6  in  the  algorithm, 
since  Lrain  and  Lmax  are  already  available  in  the  measured  function. 

3.3.3  Available  contrast 

Since  there  are  only  256  possible  inputs  d,  the  gamma  function  G(d)  is  stored  as 
a  table  of  256  luminance  values  corresponding  to  particular  digital  inputs.  We 
wish  to  arrange  a  linear  relation  between  pixel  value  p  and  screen  luminance  of 
the  form: 


L(p)  =  L0  (1  +  C  £-29)  (2) 

pO 

where  C  is  the  available  contrast,  the  maximum  contrast  that  can  be  generated 
on  a  background  L0.  Since  L(0)  >  Lmin,  we  have: 


C  s  C. 


L  -L  . 

max  mm 

L  +L  . 

mil  min 


(3) 


Cmax  is  the  maximum  possible  contrast.  If  Lmin  =  0  (i.e.  the  display  can  generate 
a  black  level  of  zero  luminance)  and  L0  is  set  to  the  mean  display  luminance 
(  =  (Lmax  +  Lmin)/2),  a  contrast  of  100%  can  be  achieved  and  Cmax  is  equal  to  1. 
In  the  linear  relation  of  eq.  (2)  pO  is  mapped  to  L0,  regardless  of  the  value  of  C. 
It  represents  a  background  level  on  which  positive  and  negative  signals  can  be 
superimposed. 


3.3.4  Inversion  of  the  gamma  function 

The  required  linear  relation  between  image  pixel  value  p  and  screen  luminance 
can  be  described  by: 

L(p)  =  L0  (1  *  Cm  (4) 


and  this  relation  is  depicted  in  the  upper  right  plane  of  Fig  7.  A  linear  mapping 
like  this  is  obtained  through  a  procedure  illustrated  in  the  same  figure.  Consider 
an  image  pixel  value  p.  By  drawing  a  straight  line  up  from  p  until  the  linear 
relation  L(p),  and  then  straight  over  to  the  ordinate  we  find  the  desired  lumi¬ 
nance  value  L(p).  In  order  to  find  which  lookup  table  input  value  yields  that 
particular  luminance,  travel  along  a  horizontal  line  at  ordinate  level  L(p)  to  the 
curve  of  the  gamma  function  G(d)  and  drop  down  to  the  abscissa.  The  point 
where  you  hit  the  d-axis  is  the  entry  that  should  be  put  into  the  correction  (look¬ 
up)  table  C  (at  address  p)  to  obtain  a  linear  output.  This  procedure  effectively 
inverts  the  gamma  function  G(d)  and  the  inverted  table  C(p)  is  also  depicted  in 
Fig.  7.  The  actual  software  procedure  is  outlined  in  the  next  section. 
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3.3.5  Procedure  for  finding  the  gumma  correction  table  C(p) 

1  Measure  G(d). 

2  For  each  p  from  1  to  255 

3  Compute  L(p)  from  Eq.  4. 

4  Find  that  d  (in  the  range  0..255)  for  which  the  value  of  G(d)  is  nearest  to 

L(p). 

5  Set  C(p)  equal  the  selected  d. 

Note  that  p  =  0  is  reserved  for  the  stimulus  background.  Therefore,  when  C(0)  is 
filled  with  [mean_red],  the  LUT  entry  for  which  the  screen  luminance  is  nearest 
to  LO,  the  gamma  correction  table  C(p)  can  the  be  used  directly  as  a  lookup 
table  to  linearize  the  display. 

3.3.6  Usable  levels 

If  the  gamma  function  were  linear,  there  would  be  a  255  distinct  luminance 
levels  (one  for  each  LUT  entry).  But  when  the  relation  is  non-linear,  the  table 
C(p)  will  have  some  levels  that  are  used  several  times,  while  others  are  skipped 
over.  This  means  that  there  are  fewer  levels  than  255  that  are  actually  used.  The 
number  of  different  levels  in  C(p)  is  called  the  number  of  usable  levels.  The 
more  nearly  linear  the  display  is,  the  higher  the  number  of  usable  levels.  The 
linearity  of  the  display  depends  on  the  settings  of  the  contrast  and  brightness 
controls  of  the  monitor.  By  carefully  adjusting  these  controls  and  after  repeated 
calibrations,  the  number  of  usable  levels  in  the  gamma  correction  table  can  be 
maximized.  To  facilitate  this  process,  the  controls  on  the  monitor  each  have 
lockable  ten-turn  potentiometers  with  a  scale.  In  practice,  we  found  that  a 
calibration  that  yields  fewer  than  170  usable  levels  (equivalent  to  about  7.4  bit) 
should  be  rejected. 


3.4  Extension  to  two  channels 

When  contrasts  lower  than  Ccrit  are  displayed,  two  channels  are  used.  In  this 
case  the  G_LUT  is  rotated  to  produce  the  required  contrast  and  the  R  LUT  is 
filled  with  a  constant  value  [red^offset]  (see  section  2.3).  The  maximum  contrast 
that  can  be  displayed  through  the  attenuated  channel  is  Ccrjt.  For  a  linear 
monitor  this  should  be  equal  to  the  attenuation  factor  for  the  green  channel  (/, 
1/16  in  our  setup),  but  for  a  non-linear  display  this  value  is  larger  by  a  factor 
equal  to  the  slope  of  the  gamma  function  in  the  point  LO.  The  exact  value  is 
determined  experimentally  as  part  of  the  calibration  procedure. 
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3.4.1  Procedure  to  determine  [ red_offsetj 

The  exact  procedure  to  determine  [red_offset]  is  as  follows: 

1  Select  a  start  value  for  [red_offset]  between  0  and  255. 

2  Load  [red_offset]  into  RJLUT. 

3  Load  the  value  of  255  into  G_LUT. 

4  Record  the  screen  luminance  as  [green_maxj. 

5  Load  the  value  of  0  into  the  G_LUT. 

6  Record  the  screen  luminance  as  [green_min], 

7  Set  [greenmean]  =  ([green_max]  +  [green_min])  /  2. 

8  If  [green  mean]  <  LO,  increase  [red_offset]  and  go  back  to  2. 

9  If  [green_mean]  >  LO,  decrease  [red_offset]  and  go  back  to  2. 

This  procedure  is  repeated  until  no  further  adjustments  in  [red_offset]  bring 

[green  mean]  any  closer  to  LO.  In  practice  we  may  assume  that  the  green 

channel  is  approximately  linear  and  therefore  we  use  a  different  procedure 
which  incorporates  this  assumption.  The  alternative  is  as  follows: 

1  Select  [meanred]  as  startvalue  for  [red^offset], 

2  Load  128  into  the  G  LUT. 

3  Repeat 

4  Load  [redoffset]  into  R  LUT. 

5  Measure  the  luminance  L. 

6  If  L  >  LO,  decrease  [red  offset], 
until  |  L  -  LO  |  is  minimal. 

3.4.2  Linearization  of  the  green  channel 

The  luminance  range  of  the  attenuated  green  channel  is  limited  to  a  small  part 
of  the  total  luminance  range  of  the  display  and  is  situated  around  LO  (when 
[red  offset]  is  loaded  into  the  R_LUT).  Tne  monitor  is  roughly  linear  in  this 
limited  range,  and  we  have  the  choice  to  either  assume  that  it  is  perfectly  linear 
or  to  calibrate  it  accurately. 

In  the  first  case  we  use  a  linear  ramp  function,  with  255  usable  levels,  as 
correction  table,  in  the  latter  case  we  measure  the  gamma  function  of  the  green 
channel  and  use  its  inversion  as  correction  table.  The  inversion  however,  is  not 
trivial  since  the  luminance  steps  in  the  attenuated  channel  are  so  small  that  they 
are  smaller  than  the  measurement  accuracy  of  our  Pritchard  photometer. 
Therefore,  the  gamma  function  of  the  green  channel  is  rather  "noisy"  and  hard  to 
invert. 

In  practice  we  assume  the  gamma  function  of  the  monitor  to  be  linear  in  the 
limited  range  and  a  255  level  ramp  function  is  used  as  correction  table  for  the 
green  channel. 
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3.4.3  Procedures  to  measure  the  critical  contrast 

There  are  two  ways  to  determine  Ccnt.  Procedure  A  is  a  simple  one  that  assumes 
that  the  screen  is  perfectly  linear  in  the  attenuated  range;  only  two  luminance 
measurements  are  used.  Procedure  B  is  more  accurate  and  takes  luminance 
measurements  over  the  whole  attenuated  range  into  account. 

Procedure  A 

1  Load  test  pattern  into  image  memory. 

2  Load  [red  offset]  into  R_LUT. 

3  Load  255  into  the  G_LUT  and  record  the  screen  luminance  as  [greenmax]. 

4  Load  0  into  the  G  LUT  and  record  the  screen  luminance  as  [green_min]. 

5  Calculate  the  critical  contrast  as 

C'crit  =  ([green_max]  -  [green_min])  /  ([green_max]  +  [green_min]) 

Procedure  B 

1  Load  test  pattern  into  image  memory. 

2  Load  [red  offset]  into  RLUT. 

3  Record  the  green  channel  Gamma  function. 

4  Fit  a  linear  regression  line  through  to  the  luminance  measurements. 

5  Calculate  the  critical  contrast  as 

[greenmin]  =  calculated  luminance  for  G_LUT  =  0 

[green  max]  =  calculated  luminance  for  GJLUT  =  255 

Ccrit  =  ([green  max]  -  [green_min])  /  ([green_max]  +  [green  min]) 

In  practice  we  use  a  modified  version  of  procedure  A  because  our  present 
Conrac  monitor  appears  to  be  almost  linear  in  the  attenuated  range.  We 
measure  the  whole  gamma  function  of  the  attenuated  channel,  not  only  the 
luminances  for  the  G  LUT  values  0  and  255,  and  extract  the  values  for 
[green  min]  and  [green_max]  from  the  data. 

3.4.4  The  number  of  usable  levels 

The  number  of  levels  N,  available  for  a  given  signal  is  proportional  to  its 
contrast  and  the  total  number  of  usable  levels  available.  For  the  two  modes  of 
operation  we  find  the  following: 

Mode  1  -  one  channel,  high  contrasts 

When  Ml  is  the  number  of  usable  levels  in  mode  1,  the  number  of  levels 
available  for  a  signal  is 
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N1  =  C  *M1 


where 


For  the  lowest  contrast  that  will  be  displayed  in  mode  1  this  comes  down  to 

Nl^nin  =  •  Ml 

In  the  present  setup  this  yields,  given  a  critical  contrast  of  about  0.1  and  190 
usable  levels,  19  distinct  levels  for  the  signal.  It  should  be  kept  in  mind  that  this 
represents  a  practical  but  worst-case  situation. 

Mode  2  -  two  channels ,  small  contrasts 

When  M2  is  the  number  of  usable  levels  in  mode  2,  the  number  of  of  levels 
available  for  a  signal  is 

N2  =  —  *  M2 
CU 


where 


0  s  C  s  C* 

cnt 


For  the  maximum  contrast  in  mode  2,  which  is  equal  to  Ccrit,  we  find:  N2  =  M2, 
which  in  the  present  setup  comes  down  to  N2  =  255.  Realize  that  in  the 
attenuated  green  channel  no  linearization  is  performed,  i.e.  no  correction  table 
with  an  inverted  Gamma  function  is  used.  For  the  smallest  contrast  of  practical 
interest  (0.005)  and  the  critical  contrast  of  our  current  setup,  the  number  of 
available  levels  still  equals  13. 

The  present  system  thus  has  19  and  13  levels  available  for  displaying  the  lowest 
possible  contrasts  in  the  two  modes,  corresponding  resp.  to  4.2  and  3.7  bits. 
These  numbers  are  low  but  are  still  considered  sufficient  for  our  purposes.  A 
possible  problem  is  that  there  is  a  large  difference  in  accuracy  between  the  lower 
end  of  mode  1  (19)  and  the  high  end  of  mode  2  (255).  A  higher  number  of 
levels  could  be  obtained  by  adding  a  third  channel,  as  was  proposed  by  Pelli  and 
Zhang  (1991).  Whether  these  figures  are  to  be  considered  sufficient  or  not 
depends  mainly  on  the  contrast  region  of  interest. 


luminance  (cd/rrv1) 
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4  TEST  RESULTS  OF  MONITOR  CALIBRATION 

In  order  to  investigate  whether  the  two  channel  contrast  control  system  performs 
accurately,  we  performed  a  test  calibration  and  with  the  results  of  that  calibra¬ 
tion  we  checked  how  wel  low  contrasts  were  rendered  on  our  display. 

The  calibration  was  done  using  our  12  bit  calibration  utility  "Calibra"  (De  Reus 
and  Valeton,  1992),  which  is  based  on  the  two  channel  contrast  control  system 
described  in  this  report.  The  gamma  function  of  our  Conrac  monitor  was 
measured  with  the  calibrated  Pritchard  photometer.  Results  are  presented  in 
Fig.  8.  As  calibration  pattern  we  used  a  512x512  pixels  uniform  field. 


63  127  191  255  63  127  191  255 


d  d 

Fig.  8  Measured  gamma  functions  for  the  (a:  left)  red  and  (b:  right) 
green  channel  in  the  two  channel  contrast  control  system. 

The  gamma  function  of  the  Conrac  monitor  is  shown  in  Fig.  8a.  The  attenuated 
green  channel  in  Fig.  8b  shows  an  almost  linear  relationship.  For  this  reason  the 
green  channel  is  assumed  to  be  linear.  The  green  gamma  function  is  rather 
"noisy"  because  the  luminance  steps  are  near  the  photometer’s  accuracy.  The 
critical  contrast  of  this  calibration  was  9.2%,  the  maximum  contrast  was  99.8%. 

To  check  the  influence  of  the  calibration  pattern  on  the  results,  w'e  did  calibra¬ 
tion  runs  with  the  test  patterns  described  in  section  3.2.1.  We  found  that  the 
calibration  pattern  had  no  significant  influence  on  the  measured  gamma  func¬ 
tions,  which  indicates  that  our  Conrac  monitor  has  a  high-voltage  power  supply 
of  excellent  quality.  Large  loads  (high  luminance  on  large  areas)  are  rendered 
almost  as  good  as  smaller  loads,  there  is  only  a  small  difference  in  maximum 
luminance. 
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The  inverted  gamma  function  of  the  red  channel  is  shown  in  Fig.  9. 


Fig.  9  Inverted  gamma  function  of  the  red  (unattenuated)  channel, 
p  is  the  image  pixel  value  and  LUT  address,  d  the  corresponding 
LUT  entry. 


desired  contrast 


Fig.  10  Desired  contrast  versus  measured  contrast  for  the  two 
channel  contrast  control  system.  Around  0.1  the  system  changes  from 
1  channel  mode  to  2  channel  mode. 


For  this  setting  of  the  monitor's  contrast  and  brightness  knobs,  the  number  of 
usable  levels  for  the  red  channel  was  190.  Since  the  green  channel  is  assumed  to 
be  linear,  the  number  of  usable  level  for  that  channel  equals  255. 


As  a  check  on  the  calibration,  the  relation  between  desired  and  displayed 
contrast  was  determined.  Contrasts  were  varied  over  3  decades  by  consecutively 
adjusting  the  pixelvalue  of  a  single  region  on  the  linearized  screen.  The  lumi¬ 
nances  of  this  region  were  measured  and  afterwards  contrasts  were  computed 
from  these  data.  The  results  are  shown  in  Fig.  10. 

From  this  figure  we  conclude  that  the  12  bit  system  works  well  since  contrasts 
over  a  broad  range  can  be  rendered  on  the  monitor  and  no  gap  between  the  two 
operation  modes  is  visible  at  the  critical  contrast  (about  0.1). 
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APPENDIX  I 


SOURCE  CODE  FOR  REAL-TIME 
CHANNEL  CONTRAST  CONTROL 


STIMULUS  DISPLAY  USING  TVVO 


SuBROo  i.NE  CMVARF21USER.  TUST.  ATTEN.  BKGINOEX,  CHAN 

>  CO.  NFRAMES.  MODE.  CRIT  CONT.  CCM1,  CCM2) 


Display  stimulus  with  time-profile  at 
M2  bit"  or  '8  o»t*  accuracy 


1  0  (original) 

2.0  (other  calling  maclumtm:  CCM1  and  CCM2) 

2.01  (changed  parameter  checking) 

April  10,  1991  -  Antoine  J.C.  da  Reus  (1.0) 

Apn)  26,  1991  -  Antoine  J.C.  do  Rous  (2.0) 

May  31.  1991  -  Antoine  J.C.  d«  Reus  (2.01) 

1  Contrary  to  me  C  version  of  this  .outme. 
the  TUST  must  be  initialized  and 
filled  with  correct  values. 

2  All  contrasts  are  translated  to  a  linear 
luminance  scale  Ov  me  arrays  CCMx  Therefore 
CCM's  are  inverted  look-up  tables  of  the 

RED  (CCM  1 )  and  GREEN  ICCM2)  channel 

3  RED  OFFSET  is  stored  in  CCM2IU 

4  This  routine  calls  the  low-level  routine 
VOCWRI1..  I  because  it  is  (much)  faster 

than  LUTRGFllO.  We  tested  VOCWRO  and  found 
that  even  3  luts  (RGB)  could  be  loaded 
m  1  field  time  (field-rate  66  Hz  = 
frame-rate  33  H*) 

6.  We  found  that  the  LUTe  are  loaded  at  the 
start  of  each  FIELD.  For  interlaced  modes, 
this  means  that  twice  ee  much  LUTe  can  be 
loaded,  than  we  expected  before  Isee  also  4.) 

0.  The  routine  is  rather  time-critical.  This  2.01 


INTEGERM  I,  J.  VOC.  VINP.  NT,  ITT.  TYNW 
INTEGERM  VOC1.  VOC2,  ITT1,  ITT2,  CHI.  CH2.  VAL 
REAL  *4  CONTRAST 

LOGICAL  MIXED  MODE 

LOGICAL  INITED.FIRST  /.TRUE./ 

REALM  TABLE<266).  RVAL 

COMMON  INITED, TABLE 

INTEGER*  2  LUTI612),  VOCBUFI32I  /32*0/.  OLDLEV,  RED  CNT 

INTEGERS  GRN  CNT,  USER2 


*  Check  variables  1 


NT  -  TLIST(I) 

IF  (NT  .LE.  0)  GOTO  1  1  1 

IF  IIATTEN  .LT.  0.0)  OR.  (ATTEN  GT  1  .0))  GOTO  1  1  1 
IF  (NFRAMES  .LE.  0)  GOTO  1  1  1 

IF  ((BKGiNDEX  .LE.  0)  OR.  (BKGINOEX  GE.  266)1  GOTO  1  1  1 
IF  ((MODE  .LT.  1)  OR  IMODE  GT.  2)1  GOTO  1  1  1 
USER2  =  USER 


If  (F(RST)  THEN 

INITED  =  FALSE. 
FIRST  -  .FALSE 
END  IF 


1  initialize  table  only  :he 
1  first  time  this  routine 
1  is  called 


IF  I.NOT  INITED)  CALL  DOINIT2 


IF  IMODE  .EQ.  1)  THEN 


c 

version  is  written  for  speed,  not  for  clsfty. 

C 

7 

Do  NOT  use  the  Fortren  function  NINTI)  to  round 

C* 

c 

integer  values,  it  will  slow  down  the  display 

c 

by  a  factor  of  two  (every  frame  instead  of 

VOC  „  0 

c 

every  field)' 

ITT  =  -1 

c 

CO  -  0 

c 

Parameters: 

VINP  =  0 

r 

USER 

11*4) 

DeAnza  user  number  10. .3) 

CALL  LUTSELIUSER,  CHAN,  CO.  ITT,  VOC) 

c 

TL1STO) 

(R*4) 

List  with  time-profile;  first 

c 

element  is  the  length  of  the  list 

C* 

c 

ATTEN 

(R  *41 

Attenuation  of  the  stimulus  (0  0-1  0) 

c 

9KG1NDEX 

(1*4) 

Background  Index 

LUTI1]  =  CCM  1 1BKGINDEXI  '  display  stimulus 

c 

CHAN 

(1*4) 

Memory  plane  m  OeArue 

OO  1  =  1. NFRAMES 

c 

CO 

11*41 

Oon’t  care 

CONTRAST  =  ATTEN  *  Tl)ST(l  ♦  1) 

c 

NFRAMES 

(1*4) 

Number  of  frames  to  display 

DO  J  =  2.266 

c 

MODE 

11*4) 

Duplay  mode: 

LUTU)  -  CcM  1  (BKGINDEX  *  NINT(CONTRAST  *  TABLEIJII) 

c 

1  ■  1  channel  (R»  "8  bit* 

END  DO 

c 

2  2  channel  iR  *  G1  "12  bit" 

CALL  VOCWRM.  USER2,  VOC8UF,  1.  1.0.  266.  LUT) 

c 

CRIT  CONT 

(R*4) 

Critical  contrast,  far  ‘12  bit’  mode 

END  DO 

c 

CCM  1 1266) 

II  *4) 

Luminance  correction  table  for  REO 

Luminance  correction  table  for  GREEN 


LUT  *  CCM  1  (BKGINOEX)  '  fdf  entire  lut  at  once 
CALL  VOCWRO,  USER?,  VOC8UF.  1.  1,0,  266.  LUT) 


IMPLICIT 

NONE 

'NTEGER’4 

USER 

1  DeAnza  user  number 

REALM 

TUSTll) 

i  time-list 

REAL  *4 

ATTEN 

1  attenuation 

INTEGERM 

BKGINOEX 

1  beckQround  index 

INTEGERM 

CHAN 

'  OeAnza  output  channel 

INTEGERM 

CO 

t  7 

INTEGERM 

NFRAMES 

1  number  of  frames 

INTEGERM 

MODE 

1  ’8  bit*  or  *12  bit* 

REALM 

CRIT  CONT 

1  criticsl  contrast 

INTEGERM 

CCM  112661 

1  correction  table  (REOl 

INTEGERM 

CCM2(266i 

'  correction  table  (GREEN) 

'  ‘12  bit*  mode 


C  Setup  video  path 

C 

C  The  VOC  section  0  >s  used  for  RED,  section  1  for  GREEN 
C  At  the  ‘12  bit’  connection,  BLUE  is  not  connected. 

C 
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CHI  - 
VOCt 
itti 
CH2  * 


CHAN 

-  0 


i 


V0C2  ^  -1 


ITT2  -  t 


•  memory  plane  for  RED  and  GREEN 
'  atart  at  VOC  Mction  0 
'  bypsM  ITT 

1  memory  plane  for  BLUE 
1  DVDM  VOC 
1  ?yp«M  ITT 


CA^L  LUTRGBSEUUSER.  CHI.  ITT1.  V0C1.  CH2.  ITT2.  VOC2) 


VOC  -  0  >UH  VOC  Motion*  0  (RED)  and  1  IGREEN) 
ITT  -  -1  i  bypass  ITT 


C*** ******  Error  Meaaagea  •*••••••••••••*••••••••••••••••• 

1  1  }  WRITEI6.  222)  ATTEN.  NFRAMES.  BKGINOEX,  NT.  MODE 
222  FORMAT C  CMVARF20:  fatal  arror:  parameter  (a)  out  of  bounds',/, 

>  I  ATTEN  ='.F6. 3.'.  NFRAMES  =  \I3.\  BKGINOEX  =  '.13. 

>  '  NT-M3.’,  MODE-MI, ’  )*) 

STOP 

C*********  Return  . . . 

333  RETURN 


. . . . . . 

c 

C  Fill  lute 
C 

C  According  to  tha  actual  contrast,  which  is  Til  ST  (I)*  ATTEN. 

C  use  one  channel  (RED)  or  two  channels  iREO  and  GREEN). 

C 

C  contrast  <  =  cm  corn  ->  GREEN.  RED  at  RED  OFFSET 
C  contrast  >  cm  corn  ->  RED,  GREEN  all  zero 


DO  I  *  1  .NFRAMES  1  frame  loop 

CONTRAST  =  ATTEN  *  TllSTlI  ♦  1 )  i  actual  contrast 

•  F  lABS(CONTRAST)  IE.  CRlT  CONT)  THEN  1  mixed  mode  or  not 
MIXED  MOOE  -  .TRUE. 

ELSE 

MIXED  MOOE  =  FALSE 
END  IF 

IF  (MIXED  MODE)  THEN 

VAL  -  CCM2I1I  RED  at  RED  OFFSET 
DO  J  -  1,266 

LUHJ>  -  VAL 
END  DO 

LUTI267)  -  CCM20  28)  <  vary  GREEN 

CONTRAST  -  CONTRAST/CRIT  CONT  i  scale-up  to  100% 
DO  J  -268.612 

VAL  -  128.0  ♦  CONTRAST  •  TABLED  -  2661  ♦  0.6 
LUTU)  -  CCM2I1  *  VAL) 

END  DO 

ELSE 

VAL  -  0  ■  GREEN  zero 

DO  J  -  267.612 
LUTU)  -  VAL 
END  DO 

LUTU)  ^  CCM 1 1BKGINOEX)  >  vary  RED 

DO  J  -  2.266 

VAL  -  BKGINOEX  *  CONTRAST  *  TABLED)  *  0.6 
LUTU)  =•  CCM  1  iVAL  *  1) 

ENO  OO 
ENO  IF 

CALL  VOCWRU.  USER2.  VOCBUF.  3.  1. 0.  612,  LUT) 

END  DO 

IF  (MIXED  MODE)  TH£N 

VAL  -  CCM20  28)  ■  reMt  GREEN  to  128 

DO  J  =  267,612 
LUTU)  -  VAL 
END  DO 

ELSE 

VAL  -  CCM  1  (BKGlNDEX)  '  reMt  RED  to  RED  OFFSET 

DO  J  -  1.266 

LUTU)  -  VAL 
END  DO 
END  IF 

CALL  VOCWRU.  USER2.  VOCBUF.  3.  1.  0,  612.  LUT) 

END  IF  i  mode  8/ 1  2  bits 


ENO 


SUBROUTINE  DOINIT2I) 


IMPLICIT  N 

INTEGER**  I 

LOGICAL  II 

REAL**  1 

COMMON  II 


NONE 


INITED 
TABLEI266) 
INITED. TABLE 


DO  I  -  0.  266 

TABLED  ♦  1)  *  I  -  128.0 
END  DO 

INITED  -  TRUE. 


RETURN 

END 


GOTO  333 
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APPENDIX  2 


CONRAC  MONITOR  SPECIFICATIONS 


This  appendix  contains  the  main  monitor  specification  of  our  CONRAC  Model 
2400  High-Resolution  Monitor.  All  figures  are  taken  from  the  " Installation , 
Operation,  and  Maintenance  Manual,  Change  1 ",  supplied  by  the  manufacturer. 


Resolution 
Geometric  distortion 

Raster  size  stability 
170  cd/m2 
Actual  display  size 
Aspect  ratio 
Video  bandwidth 
Sync 

Scan  rate  (line  frequency) 
Scan  system 


1280  horizontal  x  960  vertical  pixels 
No  point  on  raster  deviates  from  its  proper 
position  by  more  than  1%  of  raster  height 
<  1%  change  from  zero  to  100%  APL1  at 

Height:  29.46  cm,  width:  39.37  cm 
4:3  or  1:1 

100  Hz  to  40  MHz,  +  0/-3  dB 
Internal  or  external 

Any  single  scan  rate  from  15  to  37  kHz 
Interlaced  or  noninterlaced 


l 


APL  =  average  picture  level 
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